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Abstract. We present a phenomenological Lagrangian and Poisson brackets for 
obtaining nondissipative hydrodynamic theory of supersolids. A Lagrangian 
is constructed on the basis of unification of the principles of non-equilibrium 
thermodynamics and classical field theory. The Poisson brackets, governing the 
dynamics of supersolids, are uniquely determined by the invariance requirement of the 
kinematic part of the found Lagrangian. The generalization of Lagrangian is discussed 
to include the dynamics of vortices. The obtained equations of motion do not account 
for any dynamic symmetry associated with Galilean or Lorentz invariance. They can be 
reduced to the original Andreev-Lifshitz equations if to require Galilean invariance. We 
also present a relativistic-invariant supersolid hydrodynamics, which might be useful 
in astrophysical applications. 
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1. Introduction 

A supersolid state of matter was theoretically predicted by Andreev and Lifshitz in 1969 
[1]. They noticed that in crystals with a large amplitude of zero-point motion (a large 
value of the de-Boehr parameter [2]) exceptional quantum mechanical effects can occur. 
In particular, the number of lattice sites in such crystals is not equal to the number of 
atoms and vacancies exist even at the absolute zero of temperature. These vacancies are 
caused by zero-point energy, which also causes them to be mobile as waves. Since the 
vacancies are delocalized, they can be considered as weakly interacting quasiparticles. In 
a cloud of quasiparticles, there may occur a phase transition to the Bose condensed state. 
Therefore, a supersolid is a state with simultaneously broken continuous translational 
and global U{1) symmetries and, therefore, it has both crystalline and superfluid order. 
Later on, the possibility of the supersolid phase in ^He was also discussed by Chester 
[3] and Leggett [1]. For many years, in spite of numerous experimental efforts [5j, the 
supersolid behavior of ^He was not discovered. 

However, in 2004, Kim and Chan reported the possible observation of superfluidity 
in solid ^He. The superfluid density ps was observed at temperatures below about 200 
mk in a torsional oscillator containing both bulk solid helium [7j and a porous Vycor 
matrix with it [6]. The occurrence of the superfluid density in a solid helium was 
accompanied by nonclassical rotational inertia signal. Soon, several independent groups 
confirmed the evidence of nonclassical rotational inertia and the possible existence of 
the supersolid phase [8], [9], [10, [H] • At the same time the low temperature measurements 
found no evidence for superflow in solid ^He [121 [13] and some experiments [H |9] showed 
the absence of the supersolid behavior in annealed solid ^He samples. However, the 
most recent experiment showed the superfluidlike mass transport through a region with 
solid hep ^He off the melting curve Though the microscopic nature of supersolidity 
is still unclear, the experimental results on annealing confirm the recent path integral 
Monte Carlo simulations [HI [T6|, [T71 [18] that ideal crystals do not exhibit the superfluid 
properties. In addition, the simulations predict that there are no vacancies in the 
ground state of solid "^He. Note that some microscopic models for supersolids, in which 
the lattice parameters are changed independently (the lattice spacing is an independent 
thermodynamic variable), do not account for the presence of defects [T9l [20]. 

The hydrodynamic theory of supersolids was studied by a number of authors 
[U [211 ES] within the phenomenological approach and in close analogy to the two-fluid 
model. The microscopic derivation of the corresponding hydrodynamic equations was 
given in Ref. [23] on the basis of the Gibbs statistical operator and Bogolyubov's 
method of quasiaverages [2^ [25] . There exists also the useful Poisson bracket 
formalism to hydrodynamic theory of various complex systems, including liquid crystals 
[261 [271 [281 [29|, superfluid liquids [26l|30l[3ll[32l[33l[3ll[35l[36], (anti)ferromagnets 
[35] [37] . spin glasses [26], etc. Recently, an effective Lagrangian describing the low- 
energy dynamics of supersolids has been studied [38] . 

This paper concerns a phenomenological Lagrangian and the Poisson bracket 
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formalism leading to nondissipative hydrodynamics of supersolids. The main advantage 
of this formalism is that it gives the straightforward derivation of djTiamic equations 
once the dynamic variables are defined and the Poisson brackets are found. It 
also provides a tool for studying Lyapunov stability of equilibrium solutions [39] . 
Hydrodynamics of supersolids is entirely specified both by the conserved quantities 
such as the densities of momentum, entropy, particle number and the fields related to 
the broken continuous translational and global U{1) symmetries. In order to construct 
a phenomenological Lagrangian describing the low-frequency dynamics of supersolids in 
terms of the mentioned hydrodynamic fields, we have also to introduce a non-physical 
(cyclic) dynamic variable, conjugate to entropy [35]. Using the electromagnetic analogy 
with superconductors [31], in which the presence of vortices is related to the existence of 
a magnetic field, we show how to modify the obtained Lagrangian so that it reproduces 
the dynamics with vortices. In our approach, the Poisson brackets of hydrodynamic 
variables are uniquely determined by the invariance requirement of the kinematic 
part of constructed Lagrangian. The corresponding hydrodynamic equations do not 
account for any dynamic symmetry associated with Galilean or Lorentz invariance. 
The constraint on the thermodynamic potential density, following from the Galilean 
invariance, leads to Andreev-Lifshitz hydrodynamics. We also derive the relativistic- 
invariant hydrodynamics of supersolids that might be useful for describing the crystalline 
superfluidity in compact stars The limiting cases of normal and superfluid liquids 
as well as the elasticity equations are discussed. 



2. General formalism 



In this section we collect the basic results of the Hamiltonian formalism, which will 
be used below for the hydrodynamic description of superfluid systems, in particular, 
supersolids. Thus, we are interested in the slow dynamics of a set of macroscopic 
field variables T]a{x,t) = f]a{'^), where index a numbers the fields. The dynamics of 
these fields is determined by the Hamiltonian H[r)] (the Hamiltonian is a functional of 
?7a(x, t)). In order to study the evolution of macroscopic variables, consider the following 
Lagrangian: 

L = i:-H = J d'xF4x;r]]ri^i^,t)-H[r]], (1) 

where Fa[s.;ri] is a certain functional of ria{'x,t). The first term in Eq. ([T]) represents 
the kinematic part of Lagrangian. The infinitesimal transformations of the fields, 

r7a(x, t) -* ?7^(x, t) = 7]a{y^, t) + 5r7a(x, t) 

induce the variation of the action functional W = j^^ dt L, 

6W = G[t2;v]-G[h;r]] 

6H 



ti 



+ I dtj d^x'6r]fs{x',t) Q d^xJfSo.[^',x;r]]r]^{x,t) - ) , (2) 
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with 

G[t;ri] = j d^xF4x;ri]5ri^{x,t), (3) 

J.,[x,x^,]^^|>^-f4^. (4) 

The principle of stationary action {6W = 0, 5?7q,(x, ti) = ^//^(x, = 0) gives the 
following equations of motion: 

77„(x) = J rfVj-^i[x, x'; (5) 
where the inverse matrix J^^[x, x';!]] is defined by 

J rf^x" J„^[x, x"; r^] J-^i[x", x'; r/] = 6afs6{^ - x'). (6) 

Now we define the Poisson brackets for two arbitrary functionals A [77] and B[ri] as 
follows: 

where J~^[x, x';?]] = {r7Q(x), ?7/3(x')}. Then Eqs. ([5]) become 

?7„(x) = {Vai^),H} = J d'x'{riai^),r]fsi^')}j^^. (8) 

Due to the antisymmetry property Jq,/3[x, x'; r^] = — J^q,[x', x; 77], the Poisson brackets 
(I?]) meet the well-known relationships 

{A, B} = -{B, A}, {AB, C} = A{B, C} + B{A, C}, 

including the Jacobi identity 

{A {B, C}} + {B, {C, A}} + {C, {A, B}} = 0, 

which is satisfied in virtue of the following equality: 

(5Ja/3[x, x.';7]] ^ 6J^a[x.",:si;r]] ^ 5 J fs^[:>c' , x." ; rj] ^ ^ 
5?7^(x") 5'7/3(x') 5r?a(x) 

Now we address the canonical transformations. To this end, consider the 
transformations of the form: 

Va{^) =v'a{^;v], (9) 

where the transformed field rj'^i'x) is a certain functional of the initial field ria{x.). Those 
field transformations that satisfy the condition 

j d^xF^[x- r/](5r7„(x) - j d^xF^[x; v'WM = SQ[v], (10) 



or 
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are referred to as canonical transformations. Being dependent on the structure of 
canonical transformations, the quantity Q[i]] represents a functional of ?7a(x). If 
Q = const, then we deal with the homogeneous canonical transformations which leave 
invariant the kinematic part of Lagrangian (the first term in Eq. ([1])). Next, taking 
into account the identity S^Q / Sr]a{'x.)6rij3{-x') = / 5rji3{x')5rja{'^) , vahd for the second 
variational derivatives, the canonical condition f|TT]) can be written in the form 

Jc,/3 X,x;r7 = / d-'xid-'x2 ' , . t — J^a Xi , Xa ; . (12) 

J 5r?a(x) 5r//3(x') 

The Poisson brackets ([H]) are invariant under the transformations if the following 
condition is satisfied: 



J^p[^,^,V\- J d x,d X2j^^ l^i^i^^^' ^^^> 

It can be easily proved that Eq. (IT5]) is equivalent to Eq. (IT^ . which reflects the 
condition for the transformations iQ to be canonical. Indeed, introducing the notation 
Tq,/3(x, x') = 5r7^[x; ?7]/5?7/3(x'), one can write Eqs. ( |T2i) . ( |T3l) as TJ[ri']T = J[ri] and 
TJ^^[ri]T = J^^ [?]'], respectively, where Ta/3(x, x') = TaQ(x',x). Then from the first 
equation, we find J^^[ri] = T~^J''^[ri']T^^ that results in Eq. (fT3l) . ^/^"'^[r/'] = TJ~^[?7]T. 

A Lagrangian is determined up to a derivative with respect to time. It can be easily 
proved that £fc = £fc + dx{vi)/dt leads to the same Poisson brackets and equations of 
motion as £fc (see Eq. ([T])). In doing so, we have to take into account Eqs. ([5]), ( fT2l) . 

We now discuss the infinitesimal canonical transformations, 

r/„(x) r]'^{^) = r7a(x) + 6r]4:>c; r]{x.')]. (14) 

Then in case of homogeneous transformations {Q = const), the canonical condition (|TT1) 
reduces to 

/6G 
(i^x'Ja/3[x,x';r7]%[x';?7] = ^— (15) 

with 

G = J d^x'F^[^';7]]6r]f,[^';r]]. (16) 

Next, employing Eqs. ([6]), ([7]), one immediately obtains from Eq. f|T5|) : 

(5r7„[x;?7] = {?7„(x),G}. (17) 

Thus, the quantity G (see also Eq. ([3])) should be interpreted as the generator of the 
infinitesimal canonical transformations (IT^ . Such interpretation was given by Schwinger 
in the quantum action principle [HlllS]. Exactly Eqs. (|T6l) . (ITTl) will be used below to 
find the Poisson brackets of hydrodynamic variables. Finally, note that the Hamiltonian 
mechanics in terms of the arbitrary dynamic variables was studied by Pauli for the 
systems with finite degrees of freedom. 
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3. Thermodynamics 

Normal state. The normal equihbrium state of a macroscopic system can be specified 
by five independent variables: the particle number density p, the entropy density a, 
and the three components of the momentum density vTfc. Then, being a function 
of these thermodynamic variables, the energy density e determines the equation of 
state, e = e{'Kki PiO'). One can also choose another set of thermodynamic variables 
that includes the temperature T, the chemical potential /i, and the velocity Vk- 
The corresponding thermodynamic potential density is usually denoted through uj, 
10 = uj{T,fi,Vk) {uj' = ujT = —p is the Gibbs thermodynamic potential density and 
p is the pressure). 

Crystalline state. The description of the states with a spontaneously broken 
symmetry requires the introduction of the supplementary thermodynamic variables. For 
a crystal, which breaks the continuous translational symmetry, such variables specify its 
deformation. As the elastic body is deformed, each its initial coordinate i/i (Lagrangian 
coordinate) is displaced to Xj (Eulerian coordinate). These coordinates are related 
through the displacement vector Mi(x), 

Xi = yi + Ui{x.). (18) 

Differentiation of Eq. f|T8l) with respect to y^ gives 

t; — ^ii = ^ik, Ki = ^ii — t; — • (19) 
oyk OXi 

The introduced quantity Xik characterizes the crystal deformation. For an ordinary 

crystal (without defects), the number of lattice sites coincides with the number of atoms. 

Since the number of lattice sites is a constant in deformed and undeformed volumes of 

the body, n6V = no(5Vo, we come to the well-known relationship, 

n = no\\ik\, (20) 

where n, hq are, respectively, the densities of lattice sites (or atoms) in the deformed 
and undeformed volumes 6V and SVq of the body and \Xik\ denotes determinant of 
the matrix Xik. According to Eq. (120]) . the density of lattice sites n is determined 
by Xik and, therefore, it cannot be considered as the independent variable. Thus, the 
thermodynamic state of a crystal is completely specified by the three components of the 
momentum density Hk, the entropy density a, and Xik, so that e = e{7ik,(J, Xik). If the 
energy density e depends on Xik only through \Xik\, then we come back to the description 
of the normal equilibrium state (normal liquid). 

Supersolid state. As mentioned in the introduction, a supersolid state breaks both 
the continuous translational and global U{1) symmetries. Therefore, in order to describe 
the supersolid state, we have to introduce the superfluid momentum pk (along with Aj^) 
as the supplementary thermodynamic variable associated with U{1) symmetry breaking. 
The superfiuid momentum is determined by the scalar potential ip, pk = Vk'P- Besides, 
as we have already emphasized, for the supersolid phase, the density of lattice sites n {n 
is determined by Xik] see Eq. (l20l) ) does not coincide with the particle number density 
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p. Therefore, these quantities should be considered as independent thermodynamic 
variables and, consequently, e = e^Hk, p, a, Xik,Pk)- 

We now introduce another set of thermodynamic variables, more useful in 
microscopic approach to superfluid systems [231 ISSl Hlj. Moreover, as we will see 
below (see section VIII), this choice of thermodynamic variables leads to the most 
simple formulation of relativistic-invariant hydrodynamics of superfluid systems. Let us 
consider (O' = 0, k,4), p^, and Xik as new thermodynamic variables which define the 
thermodynamic potential density u = uj{Yo, Yk, Y4,pk, Xik), so that 

du = edYo + iTkdYk + pdY^ + ^r^dXik + T^dpk- (21) 

oXik dpk 

Then the densities of energy e = C,q, momentum TCk = Ck, and particle number p = (4 
are related to u by the following formulae: 

Ca = ^, a = 0,k,4. (22) 

From Eq. (!2T!) . we get the basic thermodynamic identity, 

du duj 
Yode = da - YkdiTk - Y^dp + -rr—dpk + T^dXik, 

oPk oXik 

where the entropy density is given by 

a = -uj + YaCa = -uJ + Yos + YkTTk + Y^p. (23) 
It is easy to see that the above thermodynamic identity yields 

de F4 



de 


Yo ' 


de 


Yk 




d'Kk 


Yo 


de 


1 du 


de 


1 du 


dpk 


Yo dpk ' 


dXik 


Yo dXik 



dp Yo 



(24) 



where T, f^, p are the temperature, velocity, and chemical potential, respectively. 



4. Lagrangian 



In non-equilibrium statistical mechanics, the locally equilibrium thermodynamic 
variables are assumed to describe correctly the weakly inhomogeneous states of 
condensed matter. Such a description is complete, i.e., there is no necessity to introduce 
the supplementary thermodynamic parameters in consequence of the time evolution. 

In this section, we will obtain a Lagrangian describing the supersolid dynamics 
(hydrodynamics). The starting point is the local statistical equilibrium principle. This 
principle states that on hydrodynamic time scales, the thermodynamic variables are 
assumed to be slowly varying functions of space coordinate and time. Exactly these 
functions should be considered as dynamic (hydrodynamic) variables when constructing 
the corresponding Lagrangian. In addition, the thermodynamic energy of the system 
(or the Hamiltonian) is a certain functional of the local thermodynamic variables. 

Consider the following Lagrangian density in terms of Lagrangian coordinate Ui. 

L(y) = £(y)-£(y), (25) 
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where £(y) is the kinematic part of -^^(y). The kinematic part is a hnear and 
homogeneous function in the first time derivatives of dynamic variables. In order to 
construct 'C(y), let us choose the following pairs of fields: 

7ri{y,t),Xi{y,t), a{y,t),^{y,t), p{y,t),ip{y,t) 

as conjugate dynamic variables. Some of these fields, 7ii{y,t), cr(y,t), p(y,t), are 
the local thermodynamic variables, ip{y,t) is a non-physical field representing a cyclic 
dynamic variable (the energy density e does not depend on it), (f{y,t) is a scalar field 
that defines the superfluid momentum, and Xj(y,t) is Eulerian coordinate related to 
through the displacement vector tij(x, t), 

Xi{y,t) = yi + Ui{x,t). (26) 

Note that since the formation of superfluidity breaks the global U{1) symmetry 
generated by the conserved particle number, it is natural to take the particle number 
density p(y, t) and the scalar field ip{^, t) as conjugate variables. In view of the aforesaid, 
the kinematic part density 2,{y) can be written in the form 

^{y) = My)i^iy) - ^(y)^(y) - piyMy)- (27) 

We see that £(y) is given by three terms, each of which represents the kinematic part 
associated with one or another thermodynamic degree of freedom. In particular, the 
first term is similar to the kinematic part of Lagrangian in classical mechanics, YlPiQi 
The next step is to write the Lagrangian density in terms of Eulerian coordinates 
Xi. To this end, we address Eq. ( l26l) that gives 

x,{y) = Xr,\^)u,{^) (28) 

with 

Aifc = 5ifc - VfeMi(x). (29) 

The solution Xi = Xi{y,t) of Eq. fl26|l can be inverted, yi = ?/j(x, t). For the sake 
of simplicity, we will denote Lagrangian and Eulerian fields by the same letter, e.g., 
i^iy^t) = V'(y(x?^);^) = 4'(?^:t)- Then taking into account Eq. (128!) and noting that 

^(y) = ^(x) + Xi(y)Vi?/'(x), ip{y) = ^(x) + Xi(y)Vi(/?(x), 



one obtains 



where 



SIM 



dyk 



dxi 



£(y) = g,(x)ui(x) - a(x)?/'(x) - p(x)0(x). 



gi(x) = [7rfc(x) - a(x)Vfc^/'(x) - p(x)VfeV5(x)] A^, (x 



(30) 



(31) 



and 



VTi X 



dyk 



dxi 



7ri(y), 



ax 



dyk 



dxi 



^(yj 



P(x) 



dyk 



dxi 



p(y) 



are the densities of momentum, entropy, and particle number in Eulerian coordinates, 
respectively {\dyk/dxi\ denotes Jacobian of the corresponding transformation). The 
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total Lagrangian density represents the difference between the kinematic part and the 
energy density of the system, 

L(x) = ii(x) -e [x;7ri(x'),p(x'),cr(x'),Pi(x'), Aifc(x')] . 

Here e is a certain functional of the local thermodynamic variables and Pi{x.) = Vi(f{'x) is 
the superfiuid momentum. Since the energy density must be invariant under the global 
phase transformations and spatial translations, it depends not on the quantities v'(x), 
Mj(x) but on their derivatives only, Pi(x), Ajfc(x). Note that our Lagrangian is a non- 
canonical one because ii(x) includes the nonlinear function gi(x) of dynamic variables. 
The equations of motion can be derived straightforwardly, through the principle of 
stationary action. However, we will use the Poisson-bracket approach to derive the 
dynamic equations. As we will see below, the hydrodynamic limit of these equations 
reproduces the non-dissipative Andreev-Lifshitz hydrodynamics of supersolids [1] (see 
also Refs. [211 [221 [23l [38]). 



5. Poisson brackets 



As we have already seen from the general formalism, the Poisson brackets can be 
obtained from the invariance requirement of the kinematic part of Lagrangian. Here we 
will study the transformations leaving invariant the kinematic part (!30!) and calculate 
the Poisson brackets for the set {iTi,Ui,a,ip, p,(p} of dynamic variables. We remind 
that the thermodynamic energy depends on gradients of v^(x) and Mi(x) through the 
thermodynamic quantities pj(x) = VjV9(x) and Ajfc(x) = 6ik — VfcMj(x). 

We begin with the finite time-independent transformation, Xi 
induces the following transformation law of the displacement vector: 

Ui{x) M-(X ) = Mi(x) + X- 



X X 



It 



(32) 

Equation fl32|) reflects the displacement of the physically infinitesimal volume with 
Lagrangian coordinate yi from the point Xi to x[. Since ^^(x') = 'Ut(x), 

dxi 



Aife(x) A-fc(x ) = Aii(x) 



dx'u 



(33) 



If we define the following transformation properties of our dynamic variables: 

dxi 



7r.(x) ^ vr:(x') 



o- X 



cr X 



dx'j 
dxi 



P{^) - P'(x') 



dx'j 
dxi 



dx'- 



dxk , . 



ct(x), ^(x) 
p(x), (^(x) 



^'(x) = ^(x) 



V9'(X' 



(34) 



then it is easy to show that the kinematic part fl30|) is invariant under the studied finite 
transformations. This invariance can be written in the form 



d^x (gKx)^:(x) - a'(x)^'(x) - p'(x)0'(x; 
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= j d^x (^qi{x)ui{x) - (j(x)'0(x) - p(x)0(x)j . 

Now it is evident that under the infinitesimal transformations — > = + with 
|Xi(x)| <^ 1, the variations of the dynamic variables must be given by 

5?7a(x) = <(x) - r]^{x.), r]a = {vr^, m^, a, ^p, p, (f}. 

Next, bearing in mind Eqs. fl32l) - flMl) . one obtains 

5ui(x) = Aifc(x)xfc(x), 6ip{x.) = -Xi(x)Vi^(x), (35) 

5o-(x) = -Vi(xi(x)a(x)), 6(p{x.) = -Xi(x)ViV9(x), 

5p(x) = -Vi(xi(x)p(x)), 5TTi{x.) = -Vfe(xfc(x)7ri(x)) - 7rfe(x)ViXfe(x). 

These variations preserve invariant the kinematic part 'C(x). Thus, according to the 
general formalism, they are infinitesimal canonical transformations with the following 
generator (see Eq. f|T6|) ): 



G = J d-^x (gi(x)5uj(x) - cr(x)5^/'(x) - p(x)5v9(x)) = j d-^XTTi{x.)xii^)- 

Finally, noting that the variations fl35|) assume the form 6r]a = {'705^'} (see Eq. f|T7j) ) 
and using the above explicit expression for G, one finds the Poisson brackets of 7rj(x) 
with all other dynamic variables, including the momentum density itself: 

{7rfe(x), V^(x')} = 5(x - x')VfcV^(x), {7rfe(x), or(x')} = -ct(x) Vfe(5(x - x'), 

{7rfe(x), p(x')} = -p(x) Vfe5(x - x'), {7rfe(x), (/^(x )} = 5(x - X )VA:<y9(x), 

{7ri(x), 7rA:(x')} = 7ri(x )Vfc5(x - x') - 7rfc(x)Vi5(x - x ), 

{7r,(x), Afe,(x')} = -Afe,(x) V,(5(x - x'), (36) 

where we have employed the fact that Xj(x) is an arbitrary function. 

In order to obtain the Poisson brackets missing in Eqs. ( |36l) . we have to study 
another class of the infinitesimal canonical transformations. It is easy to see that the 
field variations 

5g,(x) = 0, 5m,(x) = /,(x), 5a(x) = 0, #(x) = x(x), 5p(x) = 0, 5<^(x) = ^^(x) 

also leave invariant the kinematic part fl30l) . The arbitrary real functions /j(x), 
x(x), and ^(x) do not depend on time and dynamic variables. Following the general 
formalism, these variations can be written in terms of the Poisson brackets between the 
corresponding field and generator G (see (Eqs. [161), (|T7l)): 

{g,(x), G} = 0, K(x), G} = /,(x), {a(x), G} = 0, 

{^(x), G} = x(x), {p(x), G} = 0, {vp(x), G} = ei^), 

with 

G = [d'x (g.(x)/,(x) - a(x)x(x) - p(x)0(x)) . 



Due to the arbitrariness of the functions /j(x), x(x), and ^(x), we come to nonzero 
Poisson brackets, missing in Eqs. ( |36l) . 

{a(x'),^(x)} = 5(x-x'), {p(x'),^(x)} = 5(x-xO, (37) 
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and 



K(x),g,(x')} = M(x-^')- (38) 

Though the considered derivation gives also the vanishing Poisson brackets, we do not 
write them for the sake of brevity. Therefore, Eqs. (!36l) . (1371) provide all the nonzero 
Poisson brackets for dynamic variables {vTj, Wj, cr, p, The rest of the Poisson 
brackets between these variables turn to zero. In order to obtain all of them, it is also 
necessary to consider the kinematic part ii(x) = — gj(x)uj(x) + cr(x)^/'(x) + p(x)(y9(x), 
which differs from Eq. ( 130|) in the time derivative, and to carry out the mathematical 
manipulations similar to those described in this paragraph. Using the vanishing Poisson 
brackets and the definition of gi(x) (see Eq. (l3T|l ). it is easy to see that the bracket 
(1381) reduces to {7rj(x), Afcj(x')}, which was already found in Eqs. (|36ll . Thus, a closed 
system of the Poisson brackets is given by Eqs. fl36|) . fl37|) as well as by the vanishing 
brackets. Up to an overall minus sign, the obtained Poisson brackets agree with the 
results of Refs. [26] , where the superfluid "^He was studied. Note that the first brackets 
from Eqs. (I36l) . (1371) include a non-physical field ip{x.) introduced as a dynamic variable 
conjugate to the entropy density o"(x). We will see below that this field decouples from 
the dynamics. 

Using the Poisson brackets ( l36l) . (137|) . one can compute the bracket (j 
functionals A and B of the physical dynamic variables {vTj, Aj^, a, p, (p}, 

6A 



for 



{AB} 



SB^ SB SB^ 

op 0(f da oAkj 



5B ^ ^ SB . ^ ^ , 

OTTk 



+ 



6B 



' 6A 6A ^ 6A ^ ^ 

P^i- + cr^i^ + AfcjV 

da 



6A 



+ 



5 A 5B 

6p 6{p 



67!- i 

_ 6B6A 
6p 5ip 



6p ' 6ip UU OAkj 

where ip^i = Viip. Up to the terms with Aj^, this result coincides with the corresponding 
Poisson bracket obtained in Ref. [36] for superfluid ^He. The last term, in curly 
brackets, represents the generalized two-cocycle [36]. It comes from the Poisson bracket 
{p(x'),^(x)} = 5(x-x'). 



6. Hydrodynamic equations 

Having obtained a closed system of the Poisson brackets, we are ready to study 
the corresponding equations of motion and their long-wave limit that leads to 
supersolid hydrodynamics. Being a functional of the local thermodynamic variables, 
the Hamiltonian of the system under consideration has the form 

H = y"d3a;£[x;7r,(x'),p(x'),(T(x'),Pi(x'),Aifc(x')]. (39) 

The equations of motion can be easily obtained from Eqs. ([8]), (136|) . (1371) taking into 
account the fact that is a cyclic variable. These equations are found to be 

5H\ / 5H 6H 
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^ 6H ^ 6H ^ 5H ^ f 6H\ ^ 6H , ^ SH 

TTi + aVi- h pVi— + PiVk^ h Vfc TTi- — + TCkVi- h Afci Vj-rr— = 0, 

da op Opk V OTTfe/ OTTfc OAkj 

Ai, +V,- Aifc— =0, Pi + V, Ufc— + — =0, + ^V,V^ + — = 0, 40 

where we have used the definition of the superfiuid momentum, pi = Viip. We can see 
that ip decouples from the dynamics of the physical fields. 

We are interested in the long-wave (hydrodynamic) limit of the derived equations 
of motion. In this limit, Eqs. ( HOl) can be significantly simplified to describe 
the nondissipative hydrodynamics of supersolids. In order to obtain the desired 
hydrodynamic equations, we have to write Eqs. ( I40l) in the leading order in spatial 
gradients of the dynamic variables. It is easy to see that in zeroth order in the gradients, 
the variational derivatives of the Hamiltonian are replaced by the ordinary derivatives 
of the energy density, 

6H{r]^{^')) ^ deM^)) _ de 
5?7a(x) <9?7a(x) dr]a 

where rj^, as above, denotes the whole set of dynamic (hydrodynamic) variables. 
Therefore, in the hydrodynamic limit, Eqs. PUI) take the form 

de , de\ • „ ^, de\ 



where ji is the particle number fiux density, tik is the momentum fiux density (or the 
stress tensor), and Qi is the entropy fiux density, 

de de ^ de de ^ de de . . 

Ji = P^ + T;—^ Uk = pdik + T^i^ \-Pi^ ^Aji— — , gi = a—.{A2) 

drci dpi drck dpk dXjk drci 

The pressure p is given by 

de de de 

P=-' + ''d^+% + '''d^k ^''^ 
When obtaining the equation of motion for tTj, we have employed the constraint 
rot p = 0, which refiects the irrotational nature of the superfiuid fiow, and the evident 
property of Xik, ViXki = V^Afcj (see Eq. fl29l) ). Note that the equation for Ajfc can 
be reduced to the equation for the displacement vector pT], usually used in elasticity 
theory, 

U^-X^J^ = 0. (44) 

Equations pT]) - (l43l) provide a complete hydrodynamic description of supersolids. 
The first three equations from Eqs. (HTj) are the differential conservation laws for the 
densities of particle number, momentum, and entropy. Two other equations describe 
the time evolution of the hydrodynamic variables related to the broken symmetries (the 
field decouples from the dynamics). The conservation law for the entropy density gives 
the physical meaning of de/diii. In fact, since the superfiuid fiow is not accompanied 
by the entropy transfer, the quantity de/dni = —Yi/Yo = Vni should be interpreted as 
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the normal velocity. Eulerian coordinate Xk(y,t) (see Eq. ( 12^ ) represents the lattice 
site position of a deformed lattice and, consequently, ikiy,t) is the lattice velocity. 
Making use of Eqs. ( l28l) . (HI]), we can see that i;i(y,t) = de/diii = Vni- Therefore, the 
lattice velocity coincides with the normal velocity. Note that the obtained hydrodynamic 
equations do not account for any dynamic symmetry associated with Galilean or Lorentz 
invariance. We will study both cases below. Finally, the energy conservation law follows 
from Eqs. fj4T]) - fj43|) . Indeed, noting that 

de de de de de ■ 

£ = -T^P + TTCT + — TTi + —pi + TTT— Aifc 

dp da drvi dpi dXik 
and using Eqs.( l4Tl) -( l43l) along with the constraint rotp = 0, one obtains 

e = -VkWk, (45) 

where Wk is the energy flux density, 

de ( de de de^ de (de de^ de ^ de 



diik dp ^ da diTi ) dpk ^dp ~'~ ^* dni ) dXik dnj ' ^^^^ 

Normal liquid. Let Aj^ and Pi be the cyclic variables, so that e = e{TTi,p,a). Then 
the conservation laws for the densities of particle number, entropy, and momentum 
become (see Eqs. (iI])-(|43D, (El)) 

p + Vi{pVi) = 0, a + Vi{aVi) = 0, ni + V k{pSik + TTiVk) = 0. 

The constraint vTj = mji = mpvi, following from Galilean invariance (see below), allows 
us to transform the third equation into the Euler equation, 

Vi + {vkVk)vi = -{l/pm)Vip. 

Elastic medium. Equations of elasticity can also be obtained from Eqs. (l4ip - fH51) 
if to take into account that p and pi are the cyclic variables, i.e., e = e{TCi, ViUk, a). For 
simplicity, we restrict ourselves to the case a = 0. Then the equations for vTj and Ui 
become (we use the equation for instead of that, for Xik) 

de 



Ui = XikVk, TCi = -Vk 



(-£ + TTkVk)6ik + TTiVk + Aji- , 

dXi 



(47) 



AjfcJ 

For an ordinary crystal, not superfiuid and without defects, the density of lattice sites 
is n = no|Ajfc|. This density coincides with the particle number density p and satisfies 
the continuity equation h + Vi{nvi) = 0. Therefore, in the linear order in the deviation 
from equilibrium, the constraint TTj = mnvi gives Tij = mniii and Eq. (j47|) itself takes 
the form of Cauchy's equation, 

mnili = Vk(Tik, (48) 

where the stress tensor aik is determined by 

d'^e 

CTik = 7^ Uls = AisikUls. (49) 

duisduik 

When obtaining Eqs. (148|) . (H9l) we have employed the fact that the energy density 
depends on VjMfc through the symmetric strain tensor Uik = (l/2)(Vi'Ufc + VfcWj). The 
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elastic modulus tensor Aisik possesses the following symmetry property: Aisik = AsUk = 

^Iski ^ikls- 

Superfluid. If one consider Xik as a cyclic variable, so that e = e{nk,o; p,pk), 
then Eqs. fHT!) - (l43ll reduce immediately to equations of nondissipative superfluid 
hydrodynamics (see, e.g., Refs. [261 1^ ). 

Vortices. The developed formalism can also be modified to describe the dynamics 
of vortices in a rotating ^He. Indeed, proceeding from the analogy with superconductors 
[3T] . where the presence of vortices is related to a magnetic field, let us introduce 
a "vector potential" a. Then the local parameter uj, related to vortices, represents 
the analog of a magnetic field of superconductors, a; = V x a. In order to construct 
the corresponding Lagrangian, we also consider the "electric induction vector" d as a 
dynamic variable conjugate to a. Now the Hamiltonian H (see Eq. ( !39ll ) depends on a 
through the superfluid momentum p = Vy? — a (rot p = —a;) and d is a cyclic variable. 
The gauge- invariant kinematic part of the Lagrangian density (see Eqs. (!30l) . (!3T|) ) is 
found to be 

£(x) = gi(x)Mi(x) - cr(x)?/'(x) - p(x)(/3(x) - di(x)ai(x), 

where 

gi(x) = [7rfc(x) - o-(x)Vfc?/'(x) - p(x)(VfcV?(x) - afc(x)) - (d(x) x a;(x))fc]A^/(x). 

and the "Poynting vector" (d xw) represents the momentum density related to vortices. 
Two new dynamic variables evolve according to the following equations (we should 
remember that d is a cyclic variable): 

d = V X (v X d) — j, a = V X a;, (50) 

where j is the particle number flux density (see Eqs. ( l42i) ). Using the definition of uj 
through the "vector potential" a, one finds from the second equation 

a; = V X (v X a;). 

It can be shown that all Eqs. (l4Tl) -( l46l) . except the equation for pi, preserve their form 
in the presence of vortices [55] . The evolution of Pi is governed by 

Pi = -"^iiPkVk + /i) - (v X u)i. (51) 

Note that the relationship div d = p (the analog of Maxwell's equation) is consistent 
with the first equation from Eqs. (1501) due to the conservation law for the particle 
number density p. Thus, Eqs. ( l4T]) - (l46l) . in which the equation for pi is replaced by Eq. 
(15T1) . provide a hydrodynamic description of a superfluid (supersolid) ^He with vortices 
(see, e.g., Ref. [26]). Both quantities ip and d decouple from the dynamics. 

In conclusion of this section, we rewrite the hydrodynamic equations ( HTi) -f H6l) for 
supersolids in terms of the thermodynamic potential density u = u(YQ,Yk,Y4, Xik,Pk) 
(see section III). These equations have a very compact form and they are convenient for 
studying both Galilean and Lorentz invariance of the system. From the thermodynamic 
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relations fl22p - (l2^ . one can obtain the following form of the differential conservation 
laws for the densities of energy, momentum, and particle number: 

Ca = -VfeCafc, a = 0,z,4, (52) 

where (a = dui/dYa (Co = = C4 = P) ^"^^ Cak are the corresponding flux densities, 

Ufc - Fo dpk dYa Yo d\u dY^ Fo ^ ' 

(Cofc = Wk is the energy flux density, = Uk is the momentum flux density, and (4k = jk 
is the particle number flux density). Also, it is easy to see that the hydrodynamic 
equations for pi and Xij reduce to 

,.v.(MnP,), v,.v.(Mi), 

Equations (I52|) - (l5^ are identical to the corresponding equations of supersolid 
hydrodynamics derived within the microscopic theory [23], which is based on 
Bogolyubov's method of quasiaverages [21]. As we will show in the next section, the 
requirement of Galilean invariance leads to the Andreev-Lifshitz hydrodynamics [T]. 

7. Galilean invariance 

As we have noticed, the derived hydrodynamic equations do not account for any 
dynamic symmetry. Here we study their invariance under Galilean transformations, 
Xk — > x'l^ = Xk — Vkt, where Vk is the velocity of one inertial coordinate system with 
respect to the other. To this end, we address the hydrodynamic equations (!52|) - (!5^ in 
terms of u, since the thermodynamic potential density uj is invariant under Galilean 
transformations |44j . 

uj{Ya,Pk,Xik) = uj{Y^,p'k,X'ik)^ a = 0,A;,4, (55) 

where 

Yo ^ = Yo, Yk ^ Ffc' = Yk + VkYo, Y, ^ Y^ = Y, + mVkYk + -^Yq, 

Pk ^ p'k = Pk - mVk, Xik A-fc = Xik (56) 

and m is the mass of the ^He atom. In Eqs. (l56l) . the local thermodynamic variables 
with prime and without it are taken at the points x' and x, respectively. 

It is easy to prove the consistency of Eqs. ( l56i) with the following transformation 
law for the phase ip under Galilean transformation: 

(p ^ ip' = 1^ - mVkXk H —t. (57) 

To this end, consider the first equation from Eqs. fl54p . which assumes the form 

<P = Po, Po = ^{Y4 + YkPk). (58) 




Classical and relativistic dynamics of supersolids 



16 



Upon differentiating Eq. (1H7|) witli respect to time at fixed x' (tlie pliase </'(x) depends 
on time botli explicitly and through Xk = + V^t), we have 0' = + VfcVjtv? — 
Therefore, under Galilean transformation, po is transformed according to the law: 

Po^Po=Po + VkPk —■ 

The same transformation law is also obtained straightforwardly, by using the definition 
of Po through the hydrodynamic fields Ya and Eqs. (1561) . Thus, Eqs. (1561) . (!57j) are 
consistent. Finally, the transformation law for the displacement vector, Ui —>■ u'^ = 
Ui — Vit, leads to the invariance of Xik under Galilean transformation. 

Next, according to Eqs. fl22!) . (l23ll . the densities of energy, momentum, particle 
number, and entropy are expressed through the thermodynamic potential density u. 
Therefore, using Eqs. flS^ . fl56p . we come to the well-known transformation properties 
for these quantities under Galilean transformation, 

p p = p, a ^ a' = a 

TXk ^ Tx'k = T^k - mVkp, e e' = e - VkHk H —p- 

The transformation laws for the flux densities fl53|) have the form 

jk j'k = ik - Vkp, Uk t-fc = tik - mVijk - VkiTi + mViVkp, 

Wk w'l, = Wk- ViUk - VkS + VkViTii H ^(ifc - ^kp)- (59) 

It can be easily proved that Eqs. (!52l) - (lMl) . describing the supersolid hydrodynamics, are 
invariant under the transformations given by Eqs. (l55l) . (l56l) . ( |59ll . In this connection, 
note that the hydrodynamic fields //^(x') depend on time both explicitly and through 
the relation between the coordinates, Xk = x'^ + Vkt. 

Now let us show that the constraint on the thermodynamic potential 
density, following from Galilean invariance, leads to the Andreev-Lifshitz supersolid 
hydrodynamics. To this end, consider a reference frame K' in which the superfluid 
momentum is zero, p'^ = 0. Then Eq. fl55l) takes the form 

UJ{Ya,PkA^k)=UJ{Y:,,Q,Kk) (60) 

and the hydrodynamic fields in the reference frames K' are related, according to 
Eqs. fl56l) . by the following formulae: 

Pk = mVk = mvsk, Y^ = Yq, Y^ = Yk + VskYo, 

9 

TTIV 

Yl = Yi + mVskYk + -^^0, Kk = ^ik- 

where Vgk = Vsk{'^) is the superfluid velocity. The flux densities (155]) are expressed 
through the conserved quantities C,a = dui/dYa and dujjdpk-, dujjdXik- Using Eqs. (!60|) . 
(122]) . we can write these derivatives through their values in the reference frame K\ where 

Vsk = 0, 

2 

TTIV 

6 = e' + Vskir'k + — ^P'' TTfc = TT^ + mvskp', P = P, 

1 du ni duj duj 

+ [Vsk - VnkjP , 



Yodpk m dXik <9A- 



ik 
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Here Vnk = Vk = —Y^/Yq is the normal velocity that coincides with the lattice one. 
Therefore, the flux densities f l53|) assume the form 

Jk = \- VskP, 

m 

de 

tik = [Ta -e' + {vni - Vsi)it'i + ji' p] 5ik + w„fc7r- + v^i-n^ + mVsiVskp + ^ji^r- 



'jk 



Wk = VnkTa + VnkiVnlTr'l) + ] jk + VnlXjlT^, (61) 



where n' is the chemical potential in K', 



/i' = /i + mVniVsi (62) 

When obtaining the expressions for tj^ and Wk we have employed Eq. (123!) to eliminate 
the thermodynamic potential density uj (the entropy density is invariant under Galilean 
transformation). From the above relations, we can see the characteristic property of 
Galilean systems - the momentum density is proportional to the particle number flux 
density, mjk = i^k- Finally, remembering the thermodynamic relations ( 12^ . we can 
write Eqs. (1541) in the form 

+ + y) = 0' '^^k + ^k{\jVnj) = 0. (63) 



The derived Eqs. ( l6Tl) -( ]63l) almost coincide with the Andreev-Lifshitz equations of 
supersolid hydrodynamics [1] in the nondissipative limit. However, in contrast to 
Andreev-Lifshitz equations, the stress tensor tik includes the term nonlinear in strain 
(the last term). 

8. Lorentz invariance 

In this section we focus on the relativistic generalization of supersolid hydrodynamics. 
Instead of Galilean invariance, we will require relativistic invariance of a hydrodynamic 
theory, i.e., the invariance under the Lorentz transformations where 
= (t, x) and a'^'^ {dp^a'"^ = 6^) are the coefficients of the transformations. 
We begin the generalization of Eqs. (!52l) - (l54l) by introducing the following four- 
component quantities: 

Yf, = (Yo, Yk), = {po,Pk), >^jfi = (Ajo, Aj-fc), 

where the constraints 

Y4,+PkYk XjkYk . 

Po = 77 , Ajo = (64) 

^0 

represent the right-hand side of Eqs. (|5^ . Equations fl64l) can also be written in 
relativistic notations, 

FX = ^^4, Y.X." = 0. (65) 
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Since the volume of the system V is not a relativistic-invariant quantity, it is natural 
to introduce the relativistic-invariant potential density uj' = uj/Yq instead of uj [H] 
{uj' = —p is the Gibbs potential density, where p is the pressure) and to express the 
conserved quantities (a and their fluxes Cafc (see Eqs. ( 1221) . (153|1 ) through this relativistic- 
invariant potential density. To this end, let us consider the function a;(y^,p^, Aj^) = 
a)(Yo! '^kiPo,Pk] Ajo, Ajfc) and relate it to the Gibbs thermodynamic potential density u': 

(^'{Yo,Yk,Y4^,Pk,Xjk) = u{YQ,Yk;po,pk;XjQ,Xjk). (66) 

Next, performing the differentiation of Eq. f l66|) with respect to all hydrodynamic fields 
and taking into account the constraints one finds 

du' doj po doj Xjo duj du' duj p^ doj Xjk du 

W^^W,' Yodio ~ YodX~o' m^m^ Yodio ^ YodX~o' 
duj' dCu Yfc duj duj' 1 du) duj' dCu Y^ dCu 
dpk dpk Yodpo' dY4 Yodpo' dXjk dXjk YodXjo' 

Now, according to Eqs. fl22l) . the densities of conserved quantities can be written in 
terms of the relativistic-invariant potential density u: 

do ^ 9 duj , duj 



duj duj ^ dub 

dYk dpo dXjo 
In a similar manner, one obtains the following formulae for the corresponding flux 
densities (see Eq. flS51) ): 

_ . du _ V ^'^ doj doj 

Qk = Jk = 7^ , Cok = Wk = —^fcTTTT POT\ -^jOTTT ; 

dpk OYq dpk dXjk 

. d dtu ^ duj 

Cik = tik = -—ujYk +PiT, h Xji——. (69) 

dYi dpk dXjk 

Note that after performing the differentiation in Eqs. fl68|) . fl69l) . we should take into 
account the constraints (El]). One can easily see that the particle number density p = j° 
and the particle number flux density = jk are combined into the four-component 
vector ji" = (pjk), 

f = ^- (70) 

Similarly, the densities of energy e and momentum vr^ as well as their fluxes Wk, tik form 
the energy-momentum tensor t'^'^ of rank two, 
dujY'^ „ duj , „ duj 
oYf, dpu ^ dXj^ 

with = e, t^^ = Wk, t^^ = TTk, and t*^ = tik- Here and below, the raising and lowering 
of indices are implemented by means of the metric tensor g^^ {gQQ = 1, gik = —Sik, 
9ok = 0). 
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In terms of the introduced t^'^ and j'^, the hydrodynamic equations fl52p . (1531) assume 
the form 

d,t'^'' = 0, d,f = 0, (72) 

where = d/dx" . Equations (154|) . describing the evolution of Pi and Aj^, can also be 
written in relativistic notations, 

duP^i - dixVv = 0, d^Xiu - duXif, = 0. (73) 

Note that the first equation includes the irrotational condition of the superfiuid flow, 
rotp = 0. The second one, contains the evident property of Xik, VfcAj; = ViAjfc- Finally, 
the entropy density evolves according to the following equation: 

*-v,(.|)=o, 

which also assumes the relativistic form. Indeed, using Eqs. (1671) . one can express 
the entropy density ( l23l) through uj, a = = YoY^{duj / dY^j^) . Then and the entropy 
flux density = YkY^{duj / dY^) together form a four vector = {a^,a^) that satisfes 
the conservation law, 

d,,a'' = 0. 

So far, we have not imposed any constraints on the thermodynamic potential density 
uj associated with relativistic invariance. In fact, we have merely changed from the 
independent variables Ya (a = 0, ...4), p^, Aj^ to Y^ (/x = 0, 3), p^, Aj^, where po and 
Ajo are given by Eqs. ( l64l) . Being a relativistic invariant [H], the potential density u 
must be a function of invariants, 

UJ = uj{Ji,J2,Ji,Jij,h,Ii) (74) 

with 

Ji = (l/2)F^r^ J, = (1/2)pX, ■J^ = V,V^ J^j = (V2)A.^A/, 



Moreover, under Lorentz transformations {x^ x'^ = a^^x^), the hydrodynamic fields 
have the following transformation properties: 

Yt,^Y'^ = a^Y^, p^-^p'^ = a^p^, A^^ ^ A^-^ = a^Xj^. 

Now the energy-momentum tensor (I7T1) and the current vector fITOl) can be written in 
the final form 

and 

duj ,, duj , ,, doj ^ ,„ , 

where, after performing the necessary differentiation, we employed Eqs. (!65l) . It is easy 
to see that the energy-momentum tensor is a symmetric tensor, t^^ = f^^, as it should 
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in relativistic theory. Due to this symmetry, the momentum density coincides with the 
energy flux density. The obtained Eqs. fl72l) - fl76l) provide the hydrodynamic description 
of a relativistic supersolid. As such object, one can consider the crystalline superfluid 
state that may occur in compact stars |10]. 

Let us consider some particular cases of Eqs. (1751) . (!76ll . 
1. Relativistic normal liquid. The thermodynamic potential density u, which depends 
on Ji and ¥4^ only, describes the nondissipative hydrodynamics of a normal liquid. 



These formulae assume also the form 



where = Y^ / yY^Y^ is the four- velocity {v^v^ = 1) and w = {e + p) is the enthalpy 
density. The quantities w, e, and p are taken in the reference frame, moving with the 
normal velocity v^- 

2. Relativistic superfluid. If u does not depend on J^j and Ji, i.e., u = u)( Ji, J2, 14), 
then Eqs. (1751) . ( |76l) take the form 

The corresponding hydrodynamic equations, as well as the first equation from Eqs. (175]) . 
coincide with those obtained within the microscopic approach for a relativistic superfluid 
[44] . There is a number of equivalent formulations of nondissipative relativistic 
hydrodynamics (see, e.g., Refs. [16]- [H]) in which the authors use various hydrodynamic 
variables. However, in terms of the chosen variables, the energy-momentum tensor 
contains a smaller number of relativistic invariants. 

3. Relativistic elasticity. Let cl; be a certain function of Ji and J^j only, i.e., 
ijj = uj{Ji, Jkj,Yi). Then Eqs. fl75]) . (1751) become 

= Cog'-'^ + ^Y'^Y^ - ^A/A/, f = ^Y^. 
oJi oJkj 0Y4 



These equations along with the second equation from Eqs. (I73|) describe the relativistic 
elasticity. 

In conclusion, we have constructed a phenomenological Lagrangian that leads to 
nondissipative hydrodynamics of supersolids. We show that the Poisson brackets of 
hydrodynamic variables are found from the invariance requirement of the kinematic 
part of the constructed Lagrangian. Depending on what variables are cyclic, our 
Lagrangian and the developed approach describe not only supersolids but also normal 
and superfluid liquids, as well as the elastic bodies. We also modify our Lagrangian 
to include the dynamics of vortices in superfluid "^He. The obtained hydrodynamic 
equations of supersolids have the most general form - they do not account for Galilean 
or Lorentz invariance. The requirement of Galilean invariance gives the hydrodynamic 
equations, which almost coincide with Andreev-Lifshitz equations. The difference is that 
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our stress tensor includes the term nonlinear in strain. We also present a relativistic- 
invariant hydro dynamic theory of supersolids, which might be useful in astrophysical 
applications [IQ] . Due to the proper choice of hydrodynamic field variables, the energy- 
momentum tensor has the most simple form in comparison to other relativistic theories 
of superfluidity (it contains a smaller number of relativistic invariants; see, e.g., Refs. 
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